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>»■ Abstract 

We extract the dynamics implicit in an algebraic fitted model Hamiltonian for the hydrogen 
chromophore's vibrational motion in the molecule CF^CHFI. The original model has 4 degrees 
^ I of freedom, three positions and one representing interbond couplings. A conserved polyad al- 
■ lows the reduction to 3 degrees of freedom. For most quantum states we can identify the 
O . underlying motion that when quantized gives the said state. Most of the classifications, identi- 
2 ! fications and assignments are done by visual inspection of the already available wave function 
\ semiclassically transformed from the number representation to a representation on the reduced 
dimension toroidal configuration space corresponding to the classical action and angle variables. 



Q ' The concentration of the wave function density to lower dimensional subsets centered on ideal- 
"55 ■ ized simple lower dimensional organizing structures and the behavior of the phase along such 

^ : 

Oh- 



^ ! organizing centers already reveals the atomic motion. Extremely little computational work is 
needed. 



*present address: Center for Nonlinear and Complex Systems, Universita dell' Insubria, Via Vallegio 11, 
22100 Como, Italy 
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1 Introduction 



In recent years we have developed methods to investigate algebraic models ( spectroscopic 
Hamiltonians ) for the vibrations of molecules P3 12l El IH E] ■ The Hamiltonians reproduce and 
encode by their construction the experimental data [HI 13 IHl El HHl HH 112] • In the mentioned 
examples the system is reducible to two degrees of freedom. Therefore it is interesting to 
investigate systems where after all possible reductions 3 degrees of freedom remain in order to 
demonstrate that our methods also work well for higher dimensional systems. In this paper we 
use as example of demonstration the molecule CF3CHFI for which in Ref. |,13- an algebraic 
model has been presented for the motion of the H atom ( chromophore ). This model has 4 
degrees of freedom and it has one conserved polyad which can be used to reduce it to 3 degrees 
of freedom. 

For the already known wave function ( they are needed to construct the Hamiltonian ) trans- 
formed from the number representation so as to be able to be plotted in the toroidal configu- 
ration space of action/angle variables jHj we show that we can visually sort most of the states 
into ladders of states with similar topology. Each ladder has a relatively simple spectrum. 
Complexity arises from their interleaving. We also can recognize underlying classical lower 
dimensional organizing structures by the fact that for these inherently complex functions the 
density ( magnitude squared ) is concentrated around them and the phase has simple behavior 
near them. The position of these structures in angle space reveals the nature of the resonance 
interaction causing the topology. It also allows the reconstruction, often by analytic use of 
the transformation back to generalized coordinates, which are near displacement coordinates, 
of the motion of the atoms which underlie this particular quantum state. By counting nodes 
in plots of the wave function density and phase advances in corresponding plots of the phase 
quantum excitation numbers can be obtained. In many cases, along with the polyad quantum 
number itself, we thereby obtain a complete set of quantum numbers for a state even though 
the corresponding classical motion is nonintegrable. These classification numbers can be inter- 
preted as quasi conserved quantities for this particular state or for the ladder of states based 
on the same dynamic organizational element. Note there are usually several different types of 
organizing structures and hence ladders in the same energy region. 

Why are the wave functions that are so visually complex as to be unsortable in displacement 
coordinate space simple enough to be sorted in our semiclassical representation ? The answer is 
fivefold. First, when a spectroscopic Hamiltonian is used and as energy changes frequencies can 
move into resonance, that is rational ratio, a constant of the motion called the polyad number 
comes into existence ( see section 2 ). It can be used to reduce the dimension of the configuration 
space when the action angle variables are used, which come out in taking the semiclassical limit 
of the spectroscopic Hamiltonian. Dimension reduction is always a great simplification. Second, 
the resonant form of the fitted spectroscopic Hamiltonian can be justified formally by using 
Van Vleck or canonical perturbation theory starting with a given potential surface ^3 El 
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in a way that takes in only the dominant dynamics. The terms that are left out for the most 
part are multiple small perturbations that in turn cause multiple small distortions to the wave 
function. This makes them harder to recognize but not fundamentally different. An excellent 
example of this is seen in P, section 3 and quite pictorially in jTTj where figures 4 and 9 
show the wave functions for the same states but computed from the full and the spectroscopic 
Hamiltonian respectively. The simplicity of the latter is enough to illustrate our point. The 
third reason for the relative wave function simplicity is that as the bound wave functions lie on 
a toroidal configuration space their parts must interfere constructively as any angle is advanced 
by 27r. This coherence results in simpler but not always simple wave functions. The fourth 
reason is in the freedom to choose a particular type of reduced dimension canonical variables. 
These variables are designed with guidance from the polyad and the angle dependence of the 
resonance to have zero velocity when the phase locking caused by the resonance sets in. The 
wave function in reduced configuration space will then in turn further collapse onto the environs 
of ( in the sense of small oscillations about ) even lower dimensional surfaces or lines. We defer 
explaining and exhibiting this point till section 5 which introduces needed notation. The fifth 
reason for simplicity is that, as alluded to above, the wave function is complex and can then be 
plotted separately in density and phase. Each part helps in visually sorting of wave function 
topologies and therefore dynamics and each yields the values of different quantum numbers. 
As organizational elements we do not necessarily use actual classical periodic orbits. In many 
cases it is simpler to study the general ( average ) classical flow and to relate to it some idealized 
guiding center and to use this idealized motion as the skeletal element on which the quantum 
density is concentrated. We can imagine that this idealized structure contains the essence 
of a large number of periodic orbits performing on average a similar motion. The quantum 
dynamics is not able to resolve all the fine details of the classical motion, nor the infinity of 
long periodic orbits, nor their enormous number of bifurcations. Its limited resolution only 
recognizes the average of the classical dynamics over cells of the size of its resolution. Therefore 
we interpret the quantum state as based on the motion of the atoms which corresponds to the 
idealized structural element. In many cases one state of the group of similarly organized states 
can be interpreted as the quantum ground state on this element and the other ones are various 
longitudinal and transverse excitations of the basic structure, hence the ladders. In some cases 
only sections of a ladder are found. 

Further in this paper, we see how the higher number and smaller spacing between observed 
states of a system with more degrees of freedom makes the whole investigation more complicated 
and in particular how the higher density causes more mixing of structures. This increase of 
structural mixing might be interpreted as one way in which effects of classical chaos enter the 
quantum behavior and results in a group of states being unassignable even in the semiclassical 
angle representation that brought simplicity in appearance to the density and phases of the 
assignable states. 

In section 2 we discuss the spectroscopic Hamiltonian, section 3 presents the semiclassical 
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representation of the wave functions on the toroidal configuration space and section 4 contains 
the dimensional reduction. Section 5 explains the basic ideas how to extract the assignment of 
a wave function from its plot on the reduced configuration torus. Section 6 explains how from 
information about the dynamics in the reduced system we go back to reconstruct information 
about the original system. Section 7 contains the classification of states and section 8 gives 
final remarks. 



2 The model 

For the description of the chromophore dynamics of the molecule CF3CHFI we use the model 
set up in [T3]. It is based on 4 degrees of freedom as required by the nature of the observed 
overtone spectroscopy. Roughly they describe the three degrees of freedom of the H atom. The 
index s stands for the stretch, the index a labels the bend in the HCCF(4) plane and the index 
b labels the bend perpendicular to this plane, the index f stands for the transpositioned C-F 
stretch. For the exact atomic motion belonging to these 4 degrees of freedom see Fig. 2 in jT^. 
The algebraic Hamiltonian fitted to experiment and to calculations on a fitted potential surface 
has a natural decomposition 

Hg = Ho + W (1) 

into a diagonal part 

Ho = Yl + X! Xj^rna]ajal^a^ (2) 

(here the indices j, m run over the 4 degrees of freedom a, b, f, s) and an interaction part 

W = Y^ ksj^m[asa]aL + o.lo'jO'm]/ (2^2) + (3) 

'^[aaaaalal + a\a\aba}]/2 + (4) 
5[aaaaa}fa}f + a\a\afaf\/2 + (5) 

where aj and are the usual harmonic destruction and creation operators of degree of freedom 
j. In Eq.3 the indices j and m run over the 3 degrees of freedom a, b, f. For the coefficients we use 
the ones given in column 6 in Table I in |1^ where we see that roughly uos ~ 2uja ~ 2ujb ^ 2ujf 
giving rise to the Fermi resonances anticipated by the terms of Eq.3 and the Darling-Dennison 
resonances anticipated by the terms of Eqs.4-6. Eq.2 has first a "harmonic oscillator" term and 
a second "anharmonic" Dunham term. 

The implied unperturbed basis of the space upon which Hq operates is the number basis 

\n) = \ns,nf,na,nb) = \ns)\nf)\na)\nb) (7) 
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rij are the quantum numbers or degree of excitation of the harmonic mode j and are eigenvalues 
of the number operator fij = Ojaj. \nj) is its pure oscillator mode eigenfunction. Hq is diagonal 
on this basis. The unperturbed ( first term of i^o ) energy of a basis function |n) is E'i = J2j "^j^j- 
When frequencies are near resonant as ujs ~ 'iuja ~ 2u;b 2ujf a, polyad naturally arises with 
many consequences. First we can write E'^ = 2Pu!b with the polyad P defined as 

P = ns + {uf + na + nb)/2 (8) 

with the corresponding operator given as P by replacing nj fij in Eq.8. Clearly all basis 
functions with nj giving the same value of P are degenerate. If basis functions are sorted by 
polyad, the Hamiltonian matrix which can be evaluated algebraically, is expected by pertur- 
bation theory ( degenerate basis functions interact most strongly ) to be nearly block diagonal 
suggesting that P is a near constant of motion. In fact, it is totally block diagonal as P is here 
an exact constant of the motion since P commutes with Hg by construction albeit not with 
the true Hamiltonian to which Hg is, as said, to a high order of perturbation an approximation. 
When each block is diagonalized, each eigenfunction 

l^f ) = E 41^) (9) 

ngP 

in block P has a good conserved quantum number P which roughly measures via the equation 
for the degree of excitation of the state in units of the lowest frequency. One can now study 
each polyad separately ( second consequence of P ). Here we will concentrate on the one 161 
state polyad P = 5 which is the most complex one in the range of measured energy. 

3 Semiclassical representation 

Having the {c^} now allows us to immediately obtain, without any extra work, semiclassical 
eigenfunctions in action/angle variables ( good to the lowest two orders in h ). We do this 
because as we will see such eigenfunctions can be viewed due to the existence of the polyad as 
effectively having one less dimension. Hence the problem must now be transformed from the 
number representation to a semiclassical action/angle representation. The drop in dimension is 
characteristic of classical mechanics where a constant of the motion, here the classical analogue 
of P, can by canonical transformation allow one to reduce the dimension of the phase space. 
Here the reduction is from a four degree of freedom system with action/angle coordinates /, (p 
to a reduced system with specific value of P and three degrees of freedom with action/angle 
coordinates J, ip. This of course will accomplish the first promised simplification, that is 
dimension reduction. 

To make a transition to the corresponding classical Hamiltonian we bring the operators in Hg 
first into symmetrical order, and then apply the semiclassical substitution rules 

aj — > exp{—i(j)j),a^j I^"^ exp{i(j)j) (10) 
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where Ij and 0^ are the classical action and angle variables. The resulting classical Hamiltonian 
is using Eqs.2-6 

Ho = - l/2)ujj + ^ x,,milj - 1/2) (/^ - 1/2) (11) 

j i<m 
= E Kj,mll'^lflU^ COS(0, - 4>3 - 0m)/^ + (12) 

^IahC0s{2(t>a-2(t>,)+ (13) 

54//COs(20,-20^)+ (14) 

eV^cos(20fe-20^) (15) 

It has the conserved quantity 

K = Is+[Ia + h + If]/'2 (16) 

Because of the symmetric ordering required for the transition to classical dynamics the zero 
point between the eigenvalues of P and the corresponding value of K is shifted by 5/4. 
The reader, as an aside, might wish to contemplate how easily this transition was taken. For 
the more familiar H = T + V Hamiltonian form one must solve the Hamilton- Jacobi equation 
to get such variables; indeed a difficult process even if the system is completely integrable. The 
simplicity is due to the fact that, as said, this should be rigorously done starting with a potential 
surface and using perturbation theory. A price is paid for circumventing perturbation theory 
and using a fitted spectroscopic Hamiltonian. This price is that the relation of the a] and aj 
to the generalized momentum and position variables pj and qj is rigorously unknown. We will 
return to this issue later when going back from reduced space to displacement coordinates; our 
so called lift process. 

If one now quantizes any single harmonic oscillator term in Eqs. 11-15 using Schroedinger cor- 
respondence Ij — > —id/d(f)j, (pj — > (pj it is clear that, since the angles define a now toroidal 
configuration space, 

\nj) — > exp{inj(f)j), \n) exp{in ■ (j)) (17) 
with rij corresponding to Ij — 1/2. Prom Eq.9 we get 

*f ^ (0) = ^ 4 exp(zn • 0) (18) 

nGP 

That is, the semiclassical wave functions are already known and given by Eq.l8. 
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4 Dimension reduction 



Dimension reduction can now be carried out by noting that, using the polyad of Eq.8 to 
ehminate the fast degree of freedom s, the j-th eigenstate can be written 

*f (0) = exp(zP0,) ^P,n.,n„n, E exp(m,(0, - 0,/2)) (19) 

na,ni„nf k=a,b,f 

-exp(iP0,) J2 ^P,na,n^,nf H exp(mfc^fc) (20) 
na,nb,nf k=a,b,f 

= exp(iP0,)xf (^a, i^b, ^/) (21) 

This shows that the eigenfunctions in a given polyad P have a common phase factor dependence 
on 05 and really only depend on xf which is a function of the three angles 

ijjj = <Pj - <Ps/2 (22) 

for j = a,b, f. These three angle make up a three dimensional toroidal configuration space 
upon which Xk situated. In fact Xk J^^t the Fourier decomposition of the reduced 
dimension wave function on the configuration torus. More rigorously the wave function xf 
defined by Eqs.19- 21 actually is the semiclassical wave function after the classical constant 
of motion K is used to reduce dimension by a canonical transformation. Using the angle 
transformation of Eq.22 suggested above by the existence of P plus the trivial transformation 

e = 0. (23) 

the canonical transformation permits that the new actions are K, given by Eq.l6 and 

Jj = (24) 

again for j — a, b, f. 6 is the cyclic angle. Starting with equations 11 to 15 and using equations 
16, 22-24 a new Hamiltonian in the six canonical variables J, ip albeit parametric in P or 
can now be written so directly that to save space we omit it. The torus is the configuration 
space. Using the Schroedinger quantization for the new variables is only semiclassical, i.e. is 
correct in the two lowest orders in h. Its use gives the semiclassical expression 

for the eigenfunctions on the reduced configuration space. 

In the new coordinates the basis function with polyad number P and three indices na,ni„nf 
corresponds to the product of three plane waves or rotors as 

Yl exp(mfc^fc) (26) 

k=a,b,f 
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Each single plane wave factor loops the in the positive direction. The rotor is that into 
which the single mode harmonic oscillator has been transformed. Since Eq. 25 is a manifestly 
complex function it can be written as 

Xf W = \x^me'^^^^^ (27) 
where the phase function ^^{il^) is given as 

$f = arctan(7mxf /i^exf ) (28) 

For the basis function of Eq.26 it is well to note for future reference that the density is constant 
and the phase is 

= riftl)} + n„Va + "'fcV'b (29) 

Symmetry properties of the system in the toroidal configuration space will be used in the 
following. The Hamiltonian is invariant under a simultaneous shift of all angles by tt. 
Therefore the reduced configuration torus covers the original configuration space twice and 
all structures show up in double, even though they really exist only once. In addition the 
system is invariant under a shift of any angle ipj by 27r and it is invariant under a reflection 
where all angles go over into their negatives. 



5 Methods to evaluate wave function plots 

It is seen that the simple basis functions have constant density and any eigenstate dominated 
by a single basis function has a density without sharp localization. Resonances as for example 
2ijja ~ ojs will be seen to cause localization about a line = constant. This follows as 
■00 — 4'a ~ — constant when differentiated with respect to time gives the frequency 

relation. Hence it is seen that resonances are associated with localization and the fact that 
dipj/dt — 0. It can be claimed that by using angle coordinates that slow to zero velocity 
at resonance wc here assure that the wave function will "collapse" onto and about a lower 
dimensional subspace called the organization center. This in reverse gives a way to recognize 
the influence of resonances, namely localization of the wave function on the conflguration torus. 
Each linearly independent locking of angles therefore reduces by 1 the dimension of the subset 
of configuration space around which the wave function is concentrated. 

The fact that the wave function is on T^, a three dimensional torus, presents graphical and visual 
challenges when attempting to sort wave functions by inspection. To meet these challenges we 
employ the fact that is 1:1 with a cube with sides of angle range 27r and identifled opposite 
boundary points. That is a point on one surface of the cube is the same as a point on the direct 
opposite side of the cube. The density and phase are plotted in such cubes. The density plots 
now over the interior of the cube are organized over the whole cube or about interior planes or 



8 



lines or points if zero, one, two or three independent resonant couplings respectively are active. 
Over (or "along" in the case of a line) the organizing structure a smooth density with no sharp 
localization should appear. 

Nodal structures will be visible and countable in directions perpendicular to the organization 
element and clearly will be associated with a localized direction. The count of such transversal 
nodes supplies for each direction of localization a transverse quantum number t, that replaces 
an original mode quantum number n which has been destroyed by the resonant interaction. 
The wave function in all the localized directions can be considered as qualitatively similar 
( a continuous deformation ) of an oscillator state of the corresponding dimension. By the 
superposition of the running wave basis functions with appropriate weights and signs locally 
around the organization center the nodal pattern of an oscillator is reproduced. The transverse 
quantum numbers tk are the corresponding oscillator excitation numbers. 
The phase functions $ according to Eq.28 are smooth and close to a plane wave in the neigh- 
borhood of the organization center. They tend to have jumps by tt and singularities away from 
the center. Therefore phase advances along fundamental cycles ( loops on the torus or lines 
that connect opposite boundary points on the cube ) of the organization center are well defined 
and are necessarily integer multiples of 27r. They provide longitudinal quantum numbers 1^ 
which also replace the interacting mode quantum numbers. 

Some states present what at first look appears to be a familiar but unexpected phase pattern 
in that the phase picture is too simple. Cases will appear when the density is localized but 
the phase picture has the globally striped flag appearance of the nearly noninteracting mode 
case. The strangeness disappears if ideas from nonlinear dynamics are employed. Nonlinear 
dynamics tells us that at low P most of the motion takes place on what is called the primary 
zone of phase space. And the now interacting states are quantization of tori in this region. 
This zone has no central organizing structure other than the whole configuration space itself. 
It does contain quasiperiodic orbits corresponding to the uncoupled modes. 
As P changes coupling of a subset or even all modes can occur by at least two mechanisms. 
First is one where the existing periodic orbits bifurcate and create new secondary zones based 
on newly created organizing structures ( periodic orbits, 2 tori, 3 tori e.t.c. ). Also new 
secondary zones may be created seemingly spontaneously due to what are called saddle-center 
bifurcations. In these secondary cases the motions will have changed and states based on tori or 
near tori ( think about the resolution of h ) will have density and phase diagrams as described 
above with {P,tk,lk) quantum numbers. 

The second mechanism of coupling which will manifest itself in density localization is by con- 
tinuous deformation; a process visualized by mentally drawing the primary configuration on 
rubber sheets and stretching and pulling the sheets in various directions till the localized form 
is achieved. Under such continuous deformation the phase function $ of the continuous, single- 
valued complex wave function ^ ('(/') always must change by an integer multiple of 27r as its 
argument ip goes around any closed loop 7 in configuration space. In our case prior to defor- 
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mation ^ = n ■ tp and that integer is rij if the loop is in ipj ( this loop will be called 7^ ) i.e. as 
ipj ipj + 27r, $ ^ $ + nj27r. More generally for any $ this can be stated as the loop integral 
relation 

m,{X) = ^l^^dmX) (30) 

where rrij is an integer and A is the deformation parameter. Now under a continuous deformation 
it is clear that rrij or in our case rij must remain constant as a continuous integer valued function. 
In our case this nj can be read off from the phase diagram as equal to the number of revolutions 
of the phase along the loop jj achieved as ipj — > ipj + 27t. For the noninteracting modes the 
afore mentioned quantum number I is just the n. 

From the point of view of basis state mixing even though the phase is behaving as if a single 
basis function is totally dominant, significant state mixing can occur. This can and will occur 
in proportion to the amount of localization of the density. Weak localization tends to occur 
when one basis function is dominant, the one with the same n. 

The n assignment holds till a small change in A causes an abrupt change in that is when 
singularities occur. At this point continuous deformation will stop. Now since the wave function 
must be continuous, singularities in phase can only occur where the density and hence the wave 
function is zero. For us this occurs at the nodes and perhaps the outer edges of the localized 
oscillator. In short when nodes occur one can no longer use this simple full set of n quantum 
numbers. One can still count phase revolutions along loops inside an organized region ( parallel 
to the organizing center ) to get the / quantum numbers. Great care should be exercised in 
deciding if for a localized density the phase is a continuation of a plane wave or not. To do this 
$(^) must be a smooth function in the neighborhood of the loop along which one is counting 
the phase revolutions. 



6 Ladders and motions in displacement coordinates 

Once the wave functions in their semiclassical reduced dimension representation are diagram- 
matically presented and quantum numbers assigned as described in the previous section, as will 
be seen in the next section, most often states can be visually sorted and ordered into ladders 
of states with a given topology. Each ladder tends to lie in the Husimi sense in one of the 
phase space resonance zones discussed previously. The simplicity of sorting is due to the fact 
that as said an idealized classical organizing element underlies each ladder. This structure can 
easily be written as a configuration space relation between the V' variables. The structure will 
allow us to immediately determine what type of resonance interaction is creating this structure. 
An approximate procedure that transforms the variables ip back to displacement space reveals 
the atomic motion ( modes of motion ) whose quantization yields the ladders of states. This 
process we call "hfting" and refer to it as the "lift" . We now describe the hft in a way that 
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starts with the least approximate and most difficult to implement algorithm and proceeds to 
simpler to implement ones that still retain the qualitative spirit of the motions. 
To determine the key interactions we only illustrate the procedure with a simple example. 
Assume that the organization structure is a line ( or fiber ) along ipa = 0. Since ipa — 4>a~ 0s/2 
then 0a — 0s/2 ^ cua — ujs/^ ~ and the 2u!a — ojg phase locking relation speaks to us of 
a 2:1 Fermi resonance between the bending mode a and the stretch mode s. Other cases are 
similarly treated. 

Now for interpretation of the motion a trajectory in the (J, ■0) phase space need be found which 
when projected into configuration space should represent the organizing element. In reverse 
the knowledge of the organizing element is helpful in finding an appropriate trajectory. To 
start this search estimates of the actions will be useful. For uncoupled modes the assigned Uj 
quantum numbers using Ij = + 1/2 suffices. To obtain estimates of the initial actions the 
quantum mechanical average of the Jj can be used and trivially computed from the known 
wave functions as 

(^.) = (*|-^^l*)=E|CnlS- (31) 

The ladder state that is most localized should be used for this and to check an estimate the 
variance can be used. With this we now have initial conditions {{Jj) and a value of ipj on the 
organizing structure) with which to use Hamilton's equations to get il^j{t)i Jj{t). Then first we 
reconstruct the cyclic angle 6 by integrating 

e{t) = ^(0) + J^ds^{i^j{s), Jj{s)) (32) 

Next we undo the canonical transformation to get the corresponding (pj{t) and Ij{t). Then we 
assume some idealized harmonic model for the elementary degrees of freedom and form the 
displacement qj{t) and its conjugate momentum Pj{t) as 

Qj{t) = ^/2I~{t)cos{Mt)),pJ{t) = ^j2I~{ijsm{(l>j{t)) (33) 

The functions qj{t) will be interpreted as representatives of the atomic motion. 
Diagrams in which the qj are taken as coordinates with t as parameter give the more conven- 
tional picture of the atomic motions. In the simplified methods the average values of Ij which 
giving the scale of the qj motion will replace I{t) in Eq.33. These methods are simple because 
again Eq.33 is used but the time dependence and therefore the computation is suppressed and 
the relations gotten between the qj from the organizing structure equation that relates (pj along 
with the use of Eqs.23, 24. 

In preparation for section 7 let us discuss a few very simple cases. First there is the possibihty 

that the density is distributed rather homogeneously over the whole and the phase function 
comes close to a global plane wave with 3 component wave vector {na,nb,nf). This means no 
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locking at all, all modes run freely with their individual frequency and the excitation numbers of 
the various modes are nj read off from the wave vector for j = a,b, f (i.e. count the revolutions 
as one traverses the phase diagram in the j'th direction) and Ug = P — [ua + Ub + nf]/2. If 
the density is concentrated in a plane, e.g. the plane ipa — constant then the angle ipa does 
not move in the long run. It only may oscillate around the value of the constant. Prom the 
canonical transformation of Eq.22 we see that this means (pa — 0s/2 = constant for the angles 
of the original degrees of freedom. Therefore we have frequency and phase locking between the 
modes a and s. The other modes b and f run freely and their excitation is according to the wave 
vector of the phase function in the plane of high density. Analogous considerations hold for 
other planes of concentration of the density. If the density is concentrated along a line, e.g. the 
line ipa = constant and ipb = constant ( this is a fiber running in f direction ) then we see from 
Eq.22 that (pa — 0s/2 and (pb — 0s/2 are kept close to constant values and this means locking 
between a and s and simultaneously locking between b and s. Of course, then also a and b are 
locked relatively to each other. The mode f runs freely in this case and its excitation number 
follows from the wave number of the phase function along the fiber of high density. Analogous 
considerations will be applied for a concentration of the density along other 1 dimensional lines 
in reduced configuration space. 



7 The motion behind individual quantum states 

For polyad 5 with 161 states, the given transformations and equations lead to a three dimen- 
sional toroidal configuration space which here is a cube with coordinates ipa, i^b and ipf. Three 
dimensional ( perspective ) plots of the semiclassical state functions look too often hke large 
globes of ink and were not informative. Here we will resort to cuts of the cube ( e.g. cut 

"^a = V'fo) aiid plot density and phase in such cuts. 

The question is now where to cut. Clearly a knowledge of the classical flow was most helpful 
here. An analysis as follows is also useful. The cuts should contain organizing elements about 
which the density is locahzed. The phase should be simple over the regions of the plane having 
high density. Note the idea of an organizing structure is that one has found the configuration 
part of the quasiregular resonant zone of the coupled modes. The phase simplicity can be used 
as a tool for finding organizing structures if they exist as simple phase patterns are immediately 
apparent to the eye. Here we have used both trial and error as well as dynamics guidance to 
find the cuts. 

Dynamics guidance means for example that in CF3CHFI a resonance with a rather close fre- 
quency is Uf ujb and its associated kbj in the Hamiltonian is largest. This suggests that 
modes b and f should be coupled in the plane ipf = ipb + constant which assures 0}^-^ — ufi^^ . 
Note the constant will vary as explained below in a dynamically meaningful way, i.e. it will 
be used for state classification on a ladder of states and allow the ladder to be subdivided into 
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similar sub ladders. 

The Fermi choice uog = 2uja along with Ug = 2ujb would then give by similar reasoning the 
organizing line ipa = Ci, ipb = C2- For some states several organizing structures could be used 
and the same dynamics revealed. Our experience is to choose those corresponding resonances 
which seem more important in H and which give longer ladders of states. 

7.1 Scheme for a large number of states throughout the polyad 

Starting at the bottom of the polyad we find 64 states that definitely lie in the Darling-Dennison 
ujf = uJh resonance class in that they have densities localized about planes all denotable as 
ipf = ipb + constant. They also have simple phase plots in these planes. At the lower end a few 
of them could also be organized by the Ug = 20;/, Ug = 2u!b Fermi resonances, i.e. about a line 
in the "a" direction; but they occur when the stretch excitation is zero and this classification 
is not very physical. 

Some of the lowest states could also be assigned for all modes by quantum numbers nj and 
correspond to continuous deformations of basis functions. Although these nj are useful in the 
lift for obtaining actions these assignment misses the phase and frequency locking and the 
localization of the density. In addition we looked very carefully for primary tori in the classical 
phase space at the corresponding energy and did not find any. Therefore EBK quantization 
cannot be applied. Hence classical dynamics shows that the assignment by rij for j = a,b, f is 
unphysical. 

In what follows we shall first discuss motion in the planes ipf = ipb + const. This leads us 
to conclude that modes "a" and "s" are uncoupled and that the Darling-Dennison resonance 
couples "f " and "b" such that mode "b" , a true mode of the Hydrogen atom motion, is driven 
by the generic source mode f which represents the effect on the Hydrogens motion of the motion 
of the rest of the molecule. At this point in our analysis the quantum numbers n^, and If+b 
will first be assigned to each state as explained by the examples below. If+b represents the 
number of quanta in the DD coupling. Clearly one more quantum number is needed for a full 
assignment and this would be the transverse one representing the degree of excitation for the 
localized motion transverse to the plane ipf = ipb^ const. In most prior papers in this series this 
quantum number, t, could be counted as the number of nodes perpendicular to the organizing 
structure. In our analysis of the CDBrClF ^Hj spectra we were unable to consistently see 
these nodes clearly but could assign t as a numerical excitation index based on the distance 
of density maxima from the ipf = ipb + const organizing plane. Here the situation was not as 
opaque as nodes could be seen. What at first was confusing was that they appeared very often 
to be organized about one of two parallel planes and at higher energy about third and fourth 
planes. At low energy this observation and the fact that pairs of states with identical (n^, n^, 
lf+bi t) assignments appeared, led to the discovery of a dynamic nearly symmetric double well 
or double "valley". Two organizing planes at -i/^/ = %IJb and ■?/'/ = V'fe + ^r, each run along and 
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bisected the transverse width of one of the valleys. The ijja — was slightly deeper making 

the ipf ^ ipf + 71 invariance be only approximate. 

For some of the lower states, the energy in the "mode f - mode b" lock is low enough that 
these states are primarily locahzed in the t/jf — t/jj, or the t/jf — t/jb + n well. As energy increases 
the states are roughly speaking "above the barrier" of the "double well" with density on both 
sides and they fall into near symmetric (+) and antisymmetric (— ) pairs (with the same Ug, Ua, 
la+b, t) reflecting the approximate invariance ipb ^ '^Pb + tt. Each member of these pairs has a 
slight density preference for one of the wells. The ladders with an even start with a + state 
and ladders with an odd start with a — state. The lowest state in any ladder of constant la 
always has most of its density near the plane ipf — ipb- 

Clearly it can be said that roughly two ladders exist. Each represents states in b-f lock but with 
motions that are vr out of phase with each other. In Table 1 this classification (±) label appears. 
It is necessary for classification as the other four quantum numbers label only pairs of states. It 
is also necessary to overcome the confusion that arose in trying to count the transverse modes. 
This count now makes sense once the organizing plane is specified and it is noted that within a 
pair the t value in column 6 is common. For higher states it is sometimes becomes difficult to 
count t. The reason for this can be traced to the fact that a slice transverse to the organizing 
plane should ideally reveal wave functions that represent oscillators. The density would then 
be the highest in the two planes running parallel and on opposite sides of the organizing plane. 
Now consider that for higher states the density is localized about both planes (albeit with 
a preference for one plane). As excitation increases then the outer lobes of the transverse 
oscillator will move toward the planes in between (recall all features appear by symmetry in 
doubles) the two original planes i.e. the planes ipf = ihij±'ri/2 on which the highest density will 
now accumulate. Hence it now make sense to count transverse nodes between these now bigger 
intensities. We label these as t'. Sometimes both t and t' assignments are possible. Figure 2c 
and 2d show the t' assignment. To see it a ipf — ipa + plane that cuts the ipf — ipi, + tt/S 
plane is chosen. 

As always the transverse quantum number indicates to which extent the coupled motion goes 
out of exact phase lock. It is a measure of the width of the phase distribution. 
A feature that further complicates and hides the true interpretation of the states is mixing 
due to the accidental degeneracies of two or three states. The mixed final eigenstates that 
result are spotted by noting their near degenerate energies and their lack of almost all of the 
above discussed features. Trial and error demixing using various weights yields states with the 
above features. The notation s/s' or s/s" /s" indicates in Table I that this particular state is a 
demixing of s + s' etc. 

The main method used for assignment, at least initially, is simply to inspect various 2D cuts of 
wave functions that hopefully exhibit recognizable localizations, nodes and plane wave struc- 
tures. This "informed" inspection and slicing of the configuration cube is aided by first inspect- 
ing the lower levels on the ladders which exhibit simpler behavior. Further insight is gained 
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by running many long classical trajectories that show here that the motion, that is the flow 
of the trajectories, is mainly parallel to the ipf = ipb plane and in addition is mainly in the 
direction of the space diagonal i.e. guided by the condition ipf = ipb = ipa- This is not totally 
unexpected as ~ a;;, ~ a;/. This effect is reflected in the observation that in the organizing 
planes ipf — ipi, and ■0/ = '06 + '?'' the wave function tend to have flbers of density running along 
the diagonal. In fact a diagonal classification could have served many states as the organizing 
structure and assignment might have been made relative to it. We chose the parallel planes as 
the organizing structures as breaks in this diagonal fibration often made assignment less than 
clear. 

The "double well" observation was first made by constructing an "accessibihty diagram" . This 
is a plot of regions of configuration space, i.e. the values of ijja, 'ipb and ipf that satisfy E = 
H{Ja, Jb, Jf,'^a,4'b,'ipf', P), at each E, for given P and any compatible actions. Around 11100 
cm~l (well below the lowest quantum state of Polyad 5 and near to the classical lower end of 
Polyad 5) it was noted that only two slabs of configuration space were accessible; these are our 
wells. At the energy of the lowest quantum state all configuration space is accessible but the 
wells act as two attractive wells that cause quantum localization above them. 
To illustrate our ideas let us consider two states, state s = 7 and s = 43 at the bottom and 
the middle of the ladder respectively. Both states are organized about ipf = iph giving a DD 
ujf = Ub expectation. State 7 is a low state with its density running up the valley in the middle 
of the well as seen in Figures Ic and la. Recall again that due to the symmetry property 
discussed at the end of the section 4 all features come in double. Figure la looks down on 
the organizing plane and figure Ic looks sideways. In Figure 2 the highly exited case discussed 
above is exhibited. The reader here can without seeing it be assured that the highest density 
is in the plane ijj f = ipb + 7t/8. Figures 2a and 2b look down on this plane of high density and 
figures 2c and 2d show a view from the side that shows transverse excitation. 
With this we proceed to extract the quantum numbers Us, ria and Ib+f for the two example 
cases of states 7 and 43. Figures 1 and 2 show that if one stays on the high density plane and 
heads say in the a direction the torus is looped during a 27r change of ipa- This means that 
the ipa niotion is free and that on the phase plot following any path with ipa ^ i'a + 27r the 
free exp(ma'0a) factor in the wave function should go to exp(ina('0a + 27r)). In Fig.lb or Fig. 2b 
respectively, as ipa ^ i^a + 27r and if we follow a line in the "a" direction one sees a phase 
advance of Att and 27r respectively. We conclude that = 2 for state 7 ( Fig.lb) and Ua = 1 
for state 43 ( Fig. 2b). Sometimes here Ua will be usefully called la to emphasize that it is a 
longitudinal quantum number obtained by following the path of highest density around the "a" 
loop. 

The phase change accompanying a loop of the b direction in figures lb and 2b by the nature 

of the cut causes a looping on the ipf = ipb + constant plane in the complete configuration 
space. This looping advances both ipf and ipb by 27r and on this plane the wave function must 
change as X) cexp(m ■ tf) — > J2cexp{i{nip + 27r(n/ + n^)) The choice of the ipf = ipb + constant 
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plane required that we sec a simple plane wave factor cxp{i{nfipf + fihiph)) which means that 
all the important basis states with large mixing coefficients must have this factor. Hence the 
wave function phase must advance as we sweep ipb in this organizing plane by 27r(n/ + nj). We 
denote this by Ib+f, I for longitudinal, i.e. motion along the organizing center as opposed to 
transverse to it. This longitudinal quantum number clearly tells us the total number of quanta 
in the Darling Dennison lock. 

For state 7 in figure lb l^j^f = 8 and for state 43 in figure 2b = 7. Using the polyad 
expression for the stretch mode we get Ug = P — {ria + lb+f)/2 = for state 7 and = 1 
for state 43. Now a partial ( three out of four needed quantum numbers ) assignment can be 
made as {P,la — na,ns) or equivalently {Ib+fja — najUs). State 7 is now (5,2,0) or (8,2,0) 
respectively and state 43 is (5, 1, 1) or (7, 1, 1) respectively as recorded in table I. In principle 
Uf and Uh no longer exist as good quantum numbers except when as explained these numbers 
represent conserved integer loop integrals as in the continuous deformation case for state 7 
where both phase diagrams ( Figs lb and Id ) look as they represent 2D plane waves. Figure 
Id then tells us that Uf, — S and Uf — 0. This is useful only for obtaining 7^ = + 1/2 = 17/2 
and If = Uf + 1/2 = 1/2 which will in turn be useful in the lift process. 

For state 43 phase simplicity only happens in the organizing plane. Fig. 2d shows no useful 
phase information. The plot shows jumps along lines and ramp singularities where the phase 
value depends on the direction of approach. 

At this point and in a similar way a full 64 states for P = 5 can be assigned by Ug and Ua- 
Actually an additional 38 states seem to resemble this picture. The resemblance is only seen 
as trends from the systematic following of the 64 other ones. For these, increasing but still 
smaller effects of many other resonances are taking their toll and they significantly degrade the 
assignable states into unassignable more heavily mixed quantum analogues of classical chaos. 
The table I indicates which states have less certain assignments. 

Subladders of common values of and (or equivalently la and If+b) can now be made. To 
complete this part of the table demixing must be used and the table indicates states which 
result from the process. Each subladder now needs classification according to (+) and (— ) and 
assignment of t and/or t'. 

By looking at the wave function amplitudes above the ipf — ipb (plane 1) and ipf — ipb + 
(plane 2) wells for each pair of states with the same {ns,nb) the (+) and (— ) assignment can 
be made. It is seen that state 1 is concentrated above plane 1 and state 2 about plane 2; the 
latter is a little higher in energy. State 3 is about plane 1. The nearly degenerate states 4 and 
5 need demixing. Then (4-5) is concentrated about plane 2, and (4+5) is about both planes. 
As described above higher excitation leads to states over both wells but with a shght favoring 
of one well. 

To see the transverse structure of the states we need plots of the density in some appropriate 
plane transverse to the plane of high density. In figure Ic we see that state 7 has t = and 
in Fig. 2c we see that state 43 has t = 1. The number f does not make sense for state 7 and 
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the number t' = 2 for state 43 is rather unclear from Fig. 2c. However the phase plot in Fig. 2d 
shows two lines of phase jumps near the lines ipa = V-'fc i and these lines of phase jumps 
indicate nodal lines of the density. Since this is a somewhat indirect indication of the value 
t' — 2, in table I the number t' — 2 for state 43 is set in brackets. 

Even at this stage where other ladders discussed below are omitted spectral complexity arises 
from the interlacing of energy levels of different subladdcrs. The roughly 100 states on the 
ladder which only demands the DD-f/b lock for its construction gives a complex spectrum 
and nontransparent energy spacings. In addition the individual ladders appear irregular by 
the mixing of states. A sequence of states which are pure excitations of a single organization 
element would form a regular ladder of energy levels looking like the spectrum of one anharmonic 
oscillator. However by the mixing the eigenstates are mixtures of various pure excitations and 
thereby their energies are shifted and thereby the ladders appear more irregular. 
The reader will note that some ladder states in table I are missing. These states could not 
be found but their disappearance is confirmed using the anharmonic models expected energy 
spacing. Formally these states he among the highly mixed ( chaotic ? ) states we discuss 
below. They are mixings of more than two dynamically identifiable states and do not show the 
near degeneracy that our demixing process identifies. Parenthetically we note that the ultimate 
demixing to dynamically recognizable states is the known inverse transform of the eigenstates 
back to the dynamically uncoupled basis functions. The magnitude of the eigenvector trans- 
formation elements would then have to be used for state classification and as said at the very 
beginning, this is impossible to do for the great majority of states. 

The actions Ij can be gotten for the modes s and a using Ij = rij -\- 1/2, in addition we have 
If + lb = li)+f + 1. Then lb = Jb and // = J/ is gotten from the state functions using Eq.31. To 
observe approximate averages of the motion of the hydrogen atom the Ij can be used in Eq.32 
for 9 to enable us to get I{t) and and then using Eq. 33 to draw qualitatively pictures 
for the various pairs {qj,pj). Actually a great deal can be said without the integration for 9. 
In {qs,(la/b/f) planes a roughly rectangular region is seen whose sides are proportional to ^/2Tl 
and ^2Ia/b/f- The hydrogen trajectories motion interior to the rectangle will be quasiperiodic 
motion. The exact trajectory depends on the initial choices of 0^ and (pa/b/f and is really not 
important. In the {qa, qb) plane organization structure leads to ellipses whose relative ranges ( 
excursion from zero displacement ) in the {qa,qb) variables is I^/ h- The ellipses eccentricity 
depends on the organizing structures phase shift, viz. = z/^b + 0a = 06 + 5 — Q 
and TT give zero angular momentum and something approaching a straight line motion while 
5 = ±7r/2 gives maximal angular momentum and elliptical motion which becomes circular if 
Ia = h- The t value is the out of phase motion and causes the trajectory to change slowly its 
eccentricity. 
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7.2 Scheme for the upper end of the polyad 

At the upper end of the polyad we take to start states 160 and 161. A cut in V'b = ( or 

any other value of ipb ) reveals as we see in Fig. 3 for state 160 what looks like two columns of 
density localized around ipa = 7^ and ipf = tt. In comparison state 161 ( not shown ) only has 
one column. A cut in the plane ipf, = constant is called to see that the columns really exist. 
Figure 3 shows that they do for state 160 and it also holds for state 161. Clearly for both these 
states the line ipa — i^f — the organizing structure. Since it loops in ■05, mode b is 
decoupled and the phase picture for both states shows ( figure 3 for state 160, state 161 not 
shown ) that as ipb '^b + 27r along the organizing line no wave fronts are crossed and hence 
there is no phase advance. Hence = = 0. For the upper end of the polyad also classical 
dynamics shows a flow in b direction. As always the flgures show all structures in two copies 
which in reality are only one structure. 

Now the ipa = "0/ = 7'" organizing structure implies that the phase locks are 0^ = 0^/2 and 
(f)f = 0s/2. Time differentiation of these relations show that the frequency locks are of 2:1 
Fermi type with the stretch mode locking with the bend mode a and the mode f. Clearly for 
the organizing structure ipa = i'f then (pa = (j) j is a. valid phase lock and cu^ = a;/ is a DD 1:1 
frequency lock. Since all three types of terms appear in the Hamiltonian we can be assured that 
all three frequencies are in lock and that the quantum numbers na, nf and ns no longer exist. 
They must be replaced by three quantum numbers of the lock one of which can be P and the 
other two can be taken as the number of nodes seen in the ip^ = cut along the antidiagonal 
as ti = 1 and along the diagonal as t2 = 0. For state 161 since one column exists ti = ^2 = 0. 
If energy decreases from the top, then the classical flow turns into the space diagonal direction 
rather soon and accordingly only the 4 highest quantum states show ( partly after demixing ) a 
clean flber structure in b direction. The flber is always sharp in the a direction but broad in the 
f direction. This demonstrates again an emerging b-f mixing in the system. The next 7 states 
still show some tendency towards flbers in b direction ( they are no longer really clean ), but 
at the same time also show the beginning of flbers in the diagonal direction. Their assignment 
as states organized along b flbers is less clear but still they continue the trends seen in the 
4 highest states ( partly after demixing ). The corresponding assignments for the 11 highest 
states are compiled in table II. 

An idea of the internal motion of the deuterium and the CF stretch can be obtained from the 
actions, phase relations and the ti and t2 values. In all {qi, qj) planes the actions again deflne 
the maximal displacement yj^lj of each qj from its equilibrium. Clearly this range is meant 
in the sense that an average is made over the values of the two variables not represented in 
the plane. The plane {qb,qj), j = (1,1,3 will show rectangularly bound quasiperiodic motion. 
In the {qs,qa/f) plane a U shaped region is swept out. Here qa/f reaches its extreme values 
as qs reaches its maximum displacement and qg sweeps its range twice for each sweep of qa/f- 
Increasing ti and/or t2 causes the hfted trajectory to oscillate in a tube about this basic U lift. 
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Of course g/ and Qa reach their maxima ( minima) in phase as the three modes s, a and f are 
phase locked. 

7.3 Patterns in the dense region of the polyad 

In the middle of the polyad there are approximately 35 states with very complicated wave 
functions which we could not classify in any scheme. It is possible that further and closer 
analysis could reveal some systematics. These states are dispersed among the states mentioned 
in subsection 7.1. But interestingly a few simple states which do not fall into any scheme 
discussed so far are dispersed in this region. In this subsection we present a few such states. 
Figure 4 shows the structure of state 82, parts a and b give density and phase respectively in the 
plane ipf = and parts c ans d show the plane ipb = 'ipa + 0.29 * 27r. Part a demonstrates some 
concentration around ipa = ±7r/2 and at the same time around = V'6±7r/ 2 with a transverse 
excitation t — 1. Plots of the phases in various planes show that planes ijjb = '4>a + constant are 
the only planes with a simple phase structure. Figure 4d shows it in the plane of highest density. 
The phase function comes very close to a plane wave with la+b = 8 and If = 2. Together with 
the above mentioned transverse quantum number t = 1 this gives the complete assignment 
of the state. In Fig. 4 notice that also the interaction between the degrees of freedom f and s 
has an effect, it concentrates the density mainly around ipf = and ipf = tt. There are a few 
other examples ( states 93, 111, 114, 122 and 136 ) which show such a clean phase function in 
planes ■06 = T/'a ± 7r/2 hut do not show a clean phase function in any other plane. There are 
approximately 5 other states where the phase function is not as clean as in figure 4d but comes 
close to a continuous deformation of a plane wave only in a plane of this type. The motion of 
all these states is a circular motion of the H atom, two orthogonal bends with relative phase 
shift locked at ±7r/2. 

Numerical results for state 80 are shown in figure 5, part a and b show density and phase 
respectively in the plane iph = ipa + ipf + /2 which is the plane of high density and thereby 
the organizing structure. Part c and d show density and phase respectively in the transverse 
plane t/jf — 0. In the plane of high density the phase function comes close to a plane wave 
with la+b — 6 and = 5. This is the only plane where the phase function is a continuous 
deformation of a plane wave. This case is something new in so far as there is a plane with a 
simple phase function and at the same time not a single Uj has a definite value. There is a 
locking involving all modes without implying a locking in any pair of modes. A few other states 
( 75, 83, 87, 118 ) fall into the same scheme, i.e. they have simple phase functions only in a 
plane of the form ipb — 'ipa + i^f + constant. The implications for the motion of the atoms are 
as follows. In the nonreduced variables the locking condition is 0^ + 0^/2 — fpa — <Pf = ±7r/2. 
This can be interpreted as a locking of the beat frequency between a pair of degrees of freedom 
to the beat frequency of the remaining pair of degrees of freedom. We have several possibilities 
to group the 4 modes into two pairs. 
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8 Conclusions 



We have shown how our previously developed methods can also be applied to systems with 3 
degrees of freedom after reduction. For 75 percent of the states the strategy works successfully. 
We investigate the wave functions constructed on a toroidal configuration space. We investigate 
on which subsets the density is concentrated and evaluate the phase function on this subset of 
high density. This allows us to obtain a set of excitation numbers for the states; in the best 
cases we obtain a complete set of excitation numbers which provide a complete assignment of 
the state. The excitation numbers can be interpreted as quasi conserved quantities for this 
particular state. 

The success of the method is based on 3 properties of our strategy which more traditional 
approaches starting from potential surfaces do not have. First, we start from an algebraic 
Hamiltonian which naturally decomposes into an unperturbed Hq and a perturbation W. Per- 
turbation theory assures us that W contains only the dominant interactions, leaving out those 
that do not affect the qualitative organization of the states. In addition having the ideal Hq re- 
sults in the classical Hamiltonian being formulated from the start in action and angle variables 
belonging to Hq. For a potential surface it is very complicated to construct an appropriate Ho 
and the corresponding action/angle variables. Our knowledge of the correct Hq allows the con- 
tinuous deformation connection of some eigenstates of H with the corresponding eigenstates 
of Hq which was useful for many states. Second in the algebraic Hamiltonian the polyad ( 
which for the real molecule is only an approximate conserved quantity ) is an exactly conserved 
quantity and allows the reduction of the number of degrees of freedom. We always deal with 
one degree of freedom less than the number of degrees of freedom on the potential surface. 
Third, our strategy depends to a large extent on having a compact configuration space with 
nontrivial homology where the phase of the wave function is important and where the phase 
advances divided by 27r along the fundamental cycles of the configuration space provide a large 
part of the excitation numbers. In part also feature 1 is responsible for the importance of the 
phases. This is an essential difference to working in usual position space. In total we have 
found a semiclassical representation of the states which is actually much simpler to view than 
the usual position space one. 

The classical dynamics of the system shows rather little regular motion and is strongly unstable. 
Then the question arises in which way the quantum states reflect anything from the classical 
chaos. We propose the following answer: In the middle of the polyad a large number of 
states do not show any simple organization pattern. For some of those states it might be that 
we simply have not found the right representation in which to see the organizational center. 
However we doubt that this is the case for many states. For many of the complicated states 
the wave function is really without a simple structure, it seems to be a complicated interference 
between several organizational structures. Then immediately comes the idea to mix these states 
with their neighboring states to extract simple structures in the way as we did it for example 
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for states 22 and 24. For most of the complicated states this method does not lead to any 
success. This indicates that underlying simple patterns are distributed over many states so 
that demixing becomes rather hopeless and maybe also meaningless. In the spirit of reference 
[T^ this indicates that the density of states is already so high that many states fall into the 
energy interval over which states are considered almost degenerate and mix strongly. This is 
a strong difference to the systems we have investigated before which could be reduced to 2 
degrees of freedom. There the density of states was much lower such that mixing only occurred 
occasionally between 2 neighboring states. This stronger mixing gives some clue as to how the 
quantum wave functions become complicated for classically chaotic systems and for sufficiently 
high quantum excitation. 

Finally there may come up the question why we do not use some of the common statistical 
procedures ( like e.g. random matrix theory ) to treat our system when effects of classical chaos 
become evident in a quantum system. The answer is: We are not satisfied with only generic 
(universal) properties of the system or with the investigation to which extent such universal 
properties are realized by our system. We aim to investigate each individual quantum state one 
after the other and to work out its individual properties in order to arrive in the best case at a 
complete set of quantum numbers and in addition at a detailed description of the atomic motion 
underlying the individual state. Previous examples have shown that our procedure provides a 
complete classification for simpler systems which can be reduced to 2 degrees of freedom. As 
the example presented in this paper shows, this program is doable with good success for a large 
number of states even in a system with originally 4 degrees of freedom, reducible only to 3 
degrees of freedom and with a complicated classical dynamics. 
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Table 1: This table gives the classification and assignment of all states organized around planes 
i/jf = + constant. The first column gives the number of the state ordered according to energy. 
If the structure only becomes clear after demixing, then it contains in addition the number(s) 
of the state(s) with which it is demixed. The second column gives the energy in units of cm~^. 
The third and fourth columns give the two longitudinal quantum numbers la and li,+f- The fifth 
column gives the approximate symmetry of the state. The sixth column gives the transverse 
quantum number t relative to the planes ipf = iph and ipf = iph + n. The seventh column 
gives the transverse quantum number t' relative to the planes ipf = ipb±TT/2. If the transverse 
structure is washed out in the figures then the corresponding transverse quantum numbers are 
put in parenthesis. For such states where the individual quantum numbers ni, and Uf exist, 
these are given in columns 8 and 9. If these phase functions show some first indications of 
developing singularities, these quantum numbers arc put in brackets. States where the pattern 
is not so clear in total but where the general trends suggest the assignment have a * in the last 
column. 
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Table 2: This table gives the classification and assignment of all states organized around fibers 
in b direction. The first column gives the number of the state ordered according to energy. If 

the structure only becomes clear after demixing, then it contains in addition the number(s) of 
the statc(s) with which it is demixed. The second column gives the energy in units of cm~^. 
The third column gives the longitudinal quantum number /;,. The fourth and fifth column give 
the transverse quantum numbers ti and ^2- 
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Fig. 1 Plots of density ( part a ) and phase ( part b ) respectively in the plane ijjf = i/jj, 
and of density ( part c ) and phase ( part d ) respectively in the plane ipa = for 
the state 7. In parts a and c darker grey means higher density. In parts b and d 
points with positive real part of the wave function are marked by a dot and points 
with a positive imaginary part of the wave function are marked by an open circle. 
All frame boundaries run from to 27r. 

Fig. 2 Plots for state 43 in the plane ■0/ — V'b + Tr/S ( parts a and b ) and in plane 
■0/ = V'o + TJ" ( parts c and d ). Otherwise this plot is done in analogy to Fig.l. 

Fig. 3 Plots for state 160 in the plane = ( parts a and b ) and in plane ipf = —ipa ( 
parts c and d ). Otherwise this plot is done in analogy to Fig.l. 

Fig. 4 Plots for state 82 in the plane ipf = { parts a and b ) and in plane ipb = '0a+O.29* 27r 
( parts c and d ). Otherwise this plot is done in analogy to Fig.l. 

Fig. 5 Plots for state 80 in the plane ipb = '?/'a + '?/'/ + 7r/2 ( parts a and b ) and in plane 
Ipf = { parts c and d ). Otherwise this plot is done in analogy to Fig.l. 
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Figure 1: Plots of density ( part a ) and phase ( part b ) respectively in the plane ipf = ipf, 
and of density ( part c ) and phase ( part d ) respectively in the plane ipa — for the state 7. 
In parts a and c darker grey means higher density. In parts b and d points with positive real 
part of the wave function are marked by a dot and points with a positive imaginary part of the 
wave function are marked by an open circle. All frame boundaries run from to 27r. 



28 




PsLa psi_a 



Figure 2: Plots for state 43 in the plane il^f = iph+ir/S ( parts a and b ) and in plane ipf = ipa+T^ 
( parts c and d ). Otherwise this plot is done in analogy to Fig.l. 
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Figure 3: Fig. 3: Plots for state 160 in the plane ■06 = ( parts a and b ) and in plane -0/ 
( parts c and d ). Otherwise this plot is done in analogy to Fig.l. 
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Figure 4: Plots for state 82 in the plane ipf = { parts a and b ) and in plane iph = t/'a + 0.29 * 27r 
( parts c and d ). Otherwise this plot is done in analogy to Fig.l. 
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Figure 5: Fig.5: Plots for state 80 in the plane iph = 4>a + '4^f + 7'"/2 ( parts a and b ) and in 
plane ■^z = ( parts c and d ). Otherwise this plot is done in analogy to Fig.l. 
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